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Abstract 

The magnetic monopole in euclidean pure SU(2) gauge theory is investigated using a 
background field method on the lattice. With Monte Carlo methods we study the mass of 
the monopole in the full quantum theory. The monopole background under the quantum 
fiuctuations is induced by imposing fixed monopole boundary conditions on the walls of 
a finite lattice volume. By varying the gauge coupling it is possible to study monopoles 
with scales from the hadronic scale up to high energies. The results for the monopole 
mass are consistent with a conjecture we made previously in a realization of the dual 
superconductor hypothesis of confinement. 
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1 Introduction 



In Ref. [1] we proposed a realization of the dual superconductor hypothesis^ [2, 3] of con- 
finement in SU(2) gauge theory. We considered the contribution of a class of magnetic 
monopole configurations to the path integral, and showed that monopole condensation 
occurs within this class, leading to a string tension. The monopoles were of the 't Hooft- 
Polyakov (HP) [6, 7] form. Assuming that the string tension a is dominated by the 
contribution of these monopoles, a variational estimate yielded the value 2.3 for the fun- 
damental ratio y^/ Ajjg. This value for the string tension was close to the results of 
Monte Carlo computations. 

The classical magnetic monopole in SU(2) gauge theory is an analogue of the HP 
monopole in the SU(2) gauge-Higgs model, in particular the Bogomol'nyi-Prasad-Sommer- 
field (BPS) [8, 9] solution. This static classical solution is characterized by an arbitrary 
(mass) scale parameter /i and its mass is given by 

Ait a , , 

g 

The role of the Higgs field in this solution is played by the timelike component A4 of the 
vector potential in the euclidean pure gauge theory. Since this monopole corresponds to 
the BPS limit of zero Higgs self-coupling of the HP monopole, it is in fact a dyon. 

In developing our model of confinement we made two assumptions about the properties 
of the monopole in the quantum theory. We proposed that quantum fiuctuations induce 
an effective mass for the A4 component of the monopole field. This would cause screening 
of the electric field of the dyon. Consequently, the quantum monopole was assumed to be 
purely magnetic at large distances. This purely magnetic nature was not expected to be 
crucial for the confinement itself, but it facilitated an approximate reduction of the SU(2) 
theory to monopole electrodynamics. The second assumption was that the mass of the 
monopole in the quantum theory can be written in the form 

M= ,^1^1 CigUAR/^)), (2) 

with gji the running coupling in the i?-scheme defined in terms of the quark-antiquark 
potential, and the function C slowly increasing as for the HP monopole, 1 < C ^ 2. We 
used gji because it naturally characterizes the interactions of electric charges (/_r/2 and 
magnetic charges Air / g^. The A-parameter in the i?-scheme is related to the As in the 
MS^ Pauli-Villars and conventional lattice schemes by A^^ = 1.048 A-^ = 0.964 Apy = 
20.78 Ap [10-13]. 

We then considered monopole configurations corresponding to an arbitrary scale /i 
in which the monopoles were separated by some minimal distance b. The configura- 
tions were supposed to be semiclassical, approximate, solutions of the effective action 
with quantum fluctations taken into account. We carried out a variational calculation 
which was facilitated by an approximate mapping to a compact U(l) lattice gauge the- 
ory with lattice distance b. The variational curve in the b - plane was given by 
fib = [(/|j(Ap//i)/(/|j(6Ap)][1.59/C((/|j(Ap//i))]. Given the assumed slow variation of C, 



^The dual superconductor explanation of confinement has been criticised [4], especially the use of the 
so-called maximal abelian gauge and the abelian projection (for a review, see [5]). We intend to comment 
on the relevance of these criticisms for our work in a separate paper. 
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the value of fib along this curve was almost constant. With increasing along the 
curve, the coupling gj^ increases, the monopole mass decreases, and beyond a critical 
scale monopole condensation occurs, leading to confinement. The string tension was 
maximized leading to = 0.916~^ at a value of g\{hKR) 7.7, giving the earlier men- 
tioned value ^/(T = 2.2 Ar = 2.3 A-^. This value of the renormalized coupling could still 
be argued to be within the semiclassical domain, it corresponds to a one/two-loop beta 
function ratio of 0.77. For a typical guess C = 1.25 this gave fib = 1.10, so that the 
condensed monopoles were closely packed. The precise value of the monopole mass was 
not so important for the value of the string tension, C = 1.59 would give fib = 1. A 
change in C is compensated by a change in the minimal distance b. 

In the present paper we investigate the monopole numerically. We perform a kind of 
background-field calculation using lattice Monte-Carlo methods. The magnetic monopole 
background is induced by imposing appropriate boundary conditions. The computations 
reported here focus on a determination of the monopole mass as a function of its scale, 
to see if our assumption about the mass was justified. In particular, we determine the 
function C, defined by Eq. (2). Apart from the relevance for our model of confinement, 
this investigation is also interesting in its own right. We hope to address the issue of 
electric screening in the future. 

The paper is organized as follows. In Section 2 we recall the classical magnetic mono- 
pole in the continuum. We discuss monopole boundary conditions on a space-time cylinder 
and the role of gauge invariance in restricting the choice of possible monopole backgrounds. 
In Section 3 the monopole is put on a hypercubic lattice. We verify that local minima 
of the lattice action correspond to continuum monopoles with small discretization errors. 
Section 4 deals with the monopole in the quantum theory. We describe our Monte Carlo 
method for calculating the monopole mass and present the results of the simulations. 
Section 5 contains our conclusions. 



2 The classical magnetic monopole 

The 't Hooft-Polyakov monopole [6, 7] is a solution to the equations of motion of SU(2) 
gauge theory coupled to a triplet Higgs field in Minkowski space. We are interested in 
the analogue of this monopole in pure SU(2) gauge theory in euclidean space-time. 

We first recall the magnetic monopole and then discuss field configurations in a finite 
volume with monopole boundary conditions. 

2.1 The monopole solutions 

The Lagrangian of SU(2) gauge theory with a triplet Higgs field is given by 

9'L = -\iF;^Y + liD.^r - ^ (V'V" - ^'Y- (3) 

The quartic Higgs potential breaks the SU(2) symmetry spontaneously to U(l) and sets 
the (mass) scale fi of the classical solutions. A finite energy solution has to approach 

Iv^l^/i, D.^^^o, f;,^o, (4) 
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fast enough at large distances. At the classical level, A is related to the masses of the 
Higgs particle and the massive vector boson by 



2A = — Y' 1 ''^w = • (5) 



m 



2 
W 



The 't Hooft-Polyakov (HP) monopole [6, 7] is a static solution to the field equations 
given by the Ansatz^ 

.r.,^. ^ 1 — K(ur) , , 

Al{x,t) = e^xi 6 

r 

a/- .^ c ^ H{lir) 
LP {X,t) = dakXk , (7) 

r 

in the "radial" or "hedgehog" gauge Aq = 0, 0{x) oc x. The monopole has magnetic 
charge —1 with respect to the residual U(l) group of gauge rotations around In terms 
of the Ansatz (6,7), the field equations reduce to differential equations for the functions 
H and K with parameter A. An analytic solution to these equations has been found only 
for A = 0. It is known as the Bogomol'nyi-Prasad-Sommerfield (BPS) solution [8, 9] and 
given by 

cosh iir , , 

smn /ir 

A>r) = (9) 
smn /ir 

For A 7^ 0, solutions have been determined numerically. Fig. 1 shows H{x)lx and 
K{x) as a function of the dimensionless parameter x = [ir, for some values of A. For 
A — > oo the function H[x)/x approaches a step function while K[x) takes a limiting form 
close to the A = 12.5 curve shown. A core size of the monopole may be deduced from the 
shape of K which determines the deviation of the magnetic energy density from l/2r'*. In 
the BPS limit A = the core size may be defined as while for A — > oo it is smaller. 
Perhaps the inverse mass in units of 47r/(/^, i.e. (/iC)~^, is a good measure for the core 
size. The mass (energy) of the monopole is given by 

M = E = ^C{2X), (10) 
9 

where C(2A) = C(m|^/ Trv^) is a slowly increasing function, C(0) — 1 < C(2A) < 1.787 — 
C(oo) [14-16]. 

For A = the Higgs field is massless. This is refiected in the asymptotic behaviour of 
the Higgs field [16, 17] (the function K[fj,r) decays oc exp(— /ir)) 

^^^^ ~ /i - - - 0(exp(-/ir)) (A = 0), (11) 
H(ur) 

' ~ /i - 0(exp[-m/jr]) (0 < A < 2) 

r 

/i - 0(exp[-/ir]) (2 < A < oo) . (12) 



^Here r = |^|, a = 1, 2, 3 is a group index and kj= 2, 3 are space indices. 
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Figure 1: H[x)/x and K[x) for A = 0.02 [niH/mw = 0.2) and A = 12.5 [niH/mw = 5.0). 
The curves have been obtained numerically. The dashed curves show the A = BPS 
solution (8,9). 

If the Higgs mass is zero there is an attractive long range interaction between a pair of 
(anti)monopoles, equal in magnitude to the magnetic interaction [18]. In the massive 
Higgs case, the Higgs interaction decays exponentially so that only the magnetic interac- 
tion survives at large distances. 

We now turn to euclidean pure SU(2) gauge theory. The static HP monopole of the 
SU(2) gauge-Higgs model can be transferred to this theory by making the substitution 
Al — > A^, — > A4, and the resulting configuration is a solution of the field equations, 

Al{x,t) = cakixr ^ ' 



r 



A\{x,t) = S^kXk^^. (13) 

r 

Since there is no Higgs potential the classical solution is given by the BPS solution (8,9). 
For the moment we will consider H and K as arbitrary functions, however. 

The configuration (13) is a special case of the spherically symmetric Ansatz for SU(2) 
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gauge theory in euclidean space"^ [19], 



.a 1 - 4>2 , </'l/c ....... 

= CaklXl + —{dak - XaXk) + AiXaXk , 

Al = AoXa, (14) 

with = (f)i{t,r), Aa = A„(t,r). It will sometimes be convenient to consider the mono- 
pole from this perspective. In terms of this Ansatz the action becomes 

s = ldtld\^{F;^r 



— / dt dr [-{D^cfif + ^r'Flp + 

g'^ J-oo Jo Z 8 4 r"^ 



where = daAp — dpAa and Da4>i = da4>i — €ijAa(j)j. The four-dimensional Yang-Mills 
theory has been transformed into a two-dimensional abelian Higgs model in a curved space 
with metric = r'^Safs- The abelian gauge symmetry in this model comes from the U(l) 
group to which the original SU(2) symmetry is broken by the Ansatz (14): the abelian 
gauge transformation exp[zc<;(r, t)] corresponds to the Ad SU(2) gauge transformation 

0(t, x) = exp[iijj[t^ r)x ■ (t/2] . (16) 

For four-dimensional gauge fields of the form (14) we define the "electromagnetic" E 
and B fields by the U(l) components of the field strength, given by the radial direction 
in SU(2), 



Er = = XaF^^ = {doA^ - diAo)x, = Foix, 

1 1 1 -</)?- 



It is clear that configurations for which — > at spatial infinity, such as the BPS 
monopole (8,9) with = 0, </)2 = A", Aq = H/r^ Ai = 0, carry negative magnetic charge. 

The BPS solution is actually self-dual. Indeed, the self-duality equations are a direct 
translation of the Bogomol'nyi equations of the HP monopole in the gauge-Higgs model. 
The exponentials in H and K decay over a distance of the order and for r ^ the 
electromagnetic fields (17) fall off radially as l/r^. Hence, the monopole carries magnetic 
and electric charge and is in fact a dyon. Nevertheless, we will keep referring to it as 
"magnetic monopole", in view of the possibility that its electric charge is screened when 
quantum fluctuations are taken into account, as discussed in Ref. [1]. 

Finally, we note that the BPS solution with H replaced by —H is also a solution. This 
can be seen from the 'charge conjugation' invariance A„ — > —Aa^ </>i — > —</>!, 4>2 4*2 
of the action (15). This new solution has equal magnetic but opposite electric charge 
compared with the monopole, cf. Fq. (17). We shall call it the "anti-electric monopole". 



^The sign of A"^ and -F^j^, as well as the two-dimensional fields Aa, Fafj and (f)i is opposite to that 
in Ref. [19]. We will use iJ,,i>,. . . for four-dimensional space-time indices, a,/3, . . . for two-dimensional 
space-time indices {r,t) and . . . for the components of (f). In two dimensions, do = dt, di = dr- 
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2.2 Boundary conditions on a space-time cylinder 



We will consider gauge field configurations in the hyper-cylinder consisting of the spatial 
sphere of radius R times the periodic time interval [0,r], with fixed spatial boundary 
conditions. In a gauge field theory, it is sufficient to fix the components of the gauge 
potential tangential to the boundary, the normal components can be left free [20]. The 
path integral is a gauge-invariant functional of the fields in the boundary. 

We will need boundary conditions for monopole and vacuum configurations. The 
choice of boundary conditions is suggested by the values of the fields in these smooth 
configurations. For the vacuum the choice is obvious, the tangential are set to zero at 
the boundary. For monopole configurations we impose the boundary conditions 

1 

K 

jlQXa . (18) 

In terms of the fields in Fq. (14), </>i = </>2 = and Aq = [Iq while Ai remains unspecified 
at the boundary. These boundary conditions are compatible with a magnetic monopole 
of the form (13). We assume [IqR large enough to neglect exponential effects. We shall 
take /io > 0. Note that the anti-electric monopole would require [Iq — > — /io in Fq. (18). 

We have introduced a subscript on the parameter /io in the boundary conditions here, 
to distinguish it from the scale [i of the monopole that may be induced. The relation 
between /io and [i depends on the asymptotic behaviour of H . If H[fj,r) and K[fj,r) are 
given by the classical BPS solution (8,9), fi and /io are related hy fiQ = fi — 1/R. If the 
l/i?-term in the large- i? behaviour of H were absent, as in the case of the HP monopole 
with massive Higgs, Fq. (12), the correspondence would be /io = /<. 

Gauss' law implies that these boundary conditions fix the magnetic charge inside the 
volume, for smooth configurations. They do not fix the electric charge, however, and 
hence are not incompatible with electric screening. 

Various symmetries play a role here. First, the boundary conditions (18) are invariant 
under the abelian gauge transformation (16) with constant cu, 

= exp[iijjx ■ (t/2] . (19) 

This is true regardless of the value of R because = (j)2 = in the boundary. This 
symmetry affects the abelian part of the gauge in which configurations appear but has no 
effect on /io and the energy. 

Next, there are gauge transformations relating different values of /io. Consider an 
arbitrary smooth gauge field configuration A^ in the hypercylinder satisfying the boundary 
conditions (18). The gauge transformation 

0(t, x) = exp[7rm— X • a] , (20) 

where periodicity in time requires n to be even, transforms A^ into a new configuration 
A'^ with boundary conditions of the form (18) but with /io replaced by /io = /io + 2TTn/T. 
The configuration A'^ is singular at the origin x = 0, although gauge invariant quantities 
such as the energy density remain smooth as x — > 0. Cutting out an infinitesimal ball 
around the singularity would avoid the singularity without affecting the energy. Fq. (20) 
with odd n is also a symmetry of the action. This 'twisted gauge transformation' is 



k, tang 



Al 
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antiperiodic, 0(t = T, x) = — 0(t = 0,x), but periodicity of the gauge fields is preserved 
since they transform in the adjoint representation of SU(2). Hence, this symmetry (20) 
relates configurations with /io values differing by the discrete translations 

/io — > /"o = /"o + 27rn/r , n G ^ , (2f ) 

but with the same energy. We shall now discuss its implications for classical monopoles 
in the space-time cylinder. 

First, consider the idealized situation R — > oo. A BPS monopole in this infinite spatial 
volume has scale /i = /io + 1/R\r^oo = /"o and energy and action 

i? = ^, ^ = ^, (22) 
9 9 

see Fig. 2 (solid line). The symmetry (2f) shifts /io to other values /ig, with action 
displayed by the rising dotted lines. Note that this symmetry also allows anti-electric 
monopoles to show up at positive values of /io, compatible with the boundary conditions 
(18). They are represented by the descending dotted lines in Fig. 2. Consequently, for 
each value of /io in Fq. (18) there are various classical monopole or anti-electric mono- 
pole configurations satisfying the boundary conditions, each with different action. From 
this figure one can easily read off which of these configurations has lowest action, for a 
particular choice of /io. 




Figure 2: Action of BPS monopole configurations allowed by the boundary conditions (18) 
as a function of the parameter fioT^ for a cylinder of infinite spatial radius R = oo and 
time extent T. The solid curve represents monopoles with scale /io. Copies of monopoles 
and anti-electric monopoles under the symmetry (20,21) are indicated by dotted lines. 
The dashed horizontal line marking the instanton action, S'inst = Syr^/t/^, is shown for 
comparison. 
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In the realistic case of finite cylinder radius R the situation is slightly different. The 
boundary conditions (f8) are compatible with a BPS monopole of scale /i = /io + l/R 
provided exp(— /io-R) is negligible, cf. Eq. (If). Such a monopole has energy E = Awfi/g^^ 
but not all of it comes from inside the cylinder. The region outside contains electric and 
magnetic energy adding up to 

1 1 1 Air 

Eont = ^ dr inr^ {—- + —-) = . (23) 

Hence the monopole energy contained inside the cylinder is 

^™ = ~ 2^ = — ~ ■ (24j 

(/"^ g^'R g'^ 

It is amusing that the energy has this form also for a finite cylinder radius i?, so that 
corrections to (24) are exponential in [IqR. 

When [IqR becomes small, the exponential corrections which were neglected in formu- 
lating the boundary conditions become important. The induced object is then no longer a 
reasonable approximation to the BPS monopole with scale /i = /io + and the energy 
increases relative to the value (24). Consider the limiting case /io = 0, for example, where 
the boundary conditions enforce H[R) = and K{R) = 0. We determined the static 
minimal-action configuration inside the cylinder for these boundary conditions by numer- 
ical minimization of the energy. The solution is given hy Aq = Ai = = [H = 0) and 
(j)2 = K as shown in Fig. 3. Its energy inside the cylinder is purely magnetic and equals 

, 4^ 1.627 , , 

9^ R 

This solution is very well approximated by a cosine profile, for which the numerical con- 
stant in this expression is 1.631, only 0.26% higher. As can be seen from Fig. 3, the 
function K vanishes linearly at the boundary and the solution does not match smoothly 
to the outside region. For instance, the energy density jumps at the boundary. Apart 
from this beauty failure, this object might pass for a purely magnetic monopole with a 
scale of order l/R. Adding a magnetic outside contribution to the energy equal to half 
the value of Fq. (23), its total energy would be given by = [Air / g'^)[2A27 / R). 

A sketch of the monopole action as a function of /io in a finite cylinder is displayed in 
Fig. 4. The value of at /io = and at the other minima is proportional to T/i?, with 
proportionality constant given by Fq. (25). For i? — > oo these minima are pushed down 
into cones and the limit curves are the ones shown in Fig. 2. 

With Fig. 4 we can discuss the feasibility of a Monte-Carlo simulation of the magnetic 
monopole. For weak coupling, quantum configurations will fiuctuate around the absolute 
minimum of the action. For a given value of the boundary parameter /io, this minimum 
is readily determined from the lowest-lying solid curve in Fig. 4. Only if < /io < tt/T 
is this configuration a monopole of scale /i = /io + 1/-R, provided fioR is sufficiently large 
that the exponential effects can be neglected. If /io > tt/T, it is a monopole or anti- 
electric monopole with scale /i depending on the unique /ig, following from the shift (21), 
for which < |/io| < tt/T. It is also clear from the figure that for each classical monopole 
with /io just below tt/T (mod 2tt /T), there is an anti-electric monopole with parameter 
2tt/T — /io which has almost the same energy. A simulation at such a value of /io may 
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Figure 3: Solution for K[r) = (j)2{r) in the cylinder of finite radius R with /io = boundary 
conditions (18) (solid curve). The function cos(7rr/2i?) is shown for comparison (dashed 
curve). 

be influenced by both classical backgrounds. However, their energy densities coincide if 
/io = tt/T. 

Although it remains to be seen, in Section 3.2, how all these considerations work 
out on a hypercubic lattice, we draw the following tentative conclusion. The dynamical 
simulations will be restricted to monopole backgrounds with < /io < tt/T. As we shall 
consider energetical aspects only, it does not matter whether we deal with a monopole or 
an anti-electric monopole, and there is a preference for /io = tt/T. Conversely, for fixed /io 
only the temperature region 1/T > fio/ir is accessible. (Recall that in our notation 1/T 
is the temperature.) Simulations of static monopoles at zero temperature are excluded. 
This appears to be sufficient for application in our description of confinement, in which 
the monopoles were supposed to be semiclassical anyway. 

3 The classical monopole on the lattice 

In this Section we translate the considerations of Section 2 to the lattice. We describe how 
a monopole gauge field configuration is set up on a hypercubic lattice, discuss boundary 
conditions and symmetries, and check that minimization of the action leads to the classical 
monopole configurations. 
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Figure 4: Inside contribution to the action of classical (anti-electric) monopole configura- 
tions and their shifted copies as a function of fioT^ for a cylinder of finite spatial radius R 
and time extent T (solid curves). The approach to the asymptotic curves (dotted lines) 
is exponentially fast as i? — > oo. 

3.1 The monopole on a hypercubic lattice 

We use an X N4 lattice, with N and N4 the number of sites in the space and time 
directions. For the radially symmetric problem we are dealing with, a cubic lattice is not 
an obvious choice. However, the conventional hypercubic lattice has the advantage that 
it is easier to translate the Monte Carlo data into physical scales. There is a wealth of 
Monte Carlo data for physical quantities on hypercubic lattices to compare with. 

The lattice is set up as illustrated in Fig. 5. Although there are N sites in each spatial 
direction, the spatial volume is a cube of size [N — l)"^a"^, where a is the lattice spacing. 
We emphasize that there is no periodicity in the spatial directions and there are no links 
sticking out of the walls of the cubic box. We shall take N to be even, so that a monopole 
in the middle of the box is located at the centre of an elementary lattice cube a^. The 
time direction is periodic and has length T = N^a. 

We begin with a discussion of the classical monopole on this lattice. The relation 
between a continuum gauge field configuration and the corresponding configuration of 
link fields on the lattice is given by 

U^{x) = V exp 

with the links in the fundamental representation, 

A, = j:Al^. (27) 

a = l 



-la 



A^[x -\- tajl) dt 



(26) 
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Figure 5: Illustration of the boundary conditions for an = 4 lattice. Links lying in 
the boundary (solid lines) are fixed, dynamic links are represented by dotted lines. For 
clarity the timelike links have been drawn somewhat shorter. 



It is important to use path-ordering (denoted by V) to ensure that the gauge symmetry 
properties of the continuum gauge field are carried over to the lattice. Two continuum 
configurations related by a gauge transformation must lead to lattice configurations that 
are exactly gauge equivalent too. It is also crucial to have the correct sign in the exponent 
of Fq. (26). This sign is related to the sign in front of the A X A term in the nonabelian 
field strength F^^ which is determined by the conventions in which the classical solution 
is written down. For the monopole of Fq. (13), with scale /i, the link matrices (26) can 
be expressed in terms of conventional integrals. 



Uk{x,t) 

Ui{x^ t) 
r{s) 



exp 



exp 



2 Jo 



1 1 - K{nr{s)) 



ds 



2 ' 



H{tir) 



(28) 
(29) 



r(0). 



In order to get an idea about discretization errors we compute the energy of the 
classical BPS solution (8,9), latticized in this way. In this context we think of the energy 
as the action per unit of time. The energy of the lattice configuration inside the box is 
calculated with the plaquette action 



^ = 4 E (l-^tr[/,.(x 



(30) 



It has to be compared with the continuum result for the energy inside the cubic box. 

To get accurate results, we have to assign relative weights < 1 to the plaquettes in the 
boundary of the lattice volume, as can be seen as follows. The lattice action is defined as 
the sum of the plaquettes, up to a proportionality constant. However, for our finite lattice 
volume a more precise definition of the total action is the sum of the action contained 
in each of the N4[N — 1)"^ elementary hypercubes lying inside the volume. The action 



12 



contained in one a"^ elementary hypercube, and assigned to the centre of this hypercube, 
is I times the sum of the 24 plaquettes framing it. The factor | comes from the fact that 
each plaquette contributes to 4 hypercubes. If a plaquette lies in the boundary of the 
lattice box, however, only 1 or 2 of its 4 hypercubes lie inside the box. Such a plaquette 
has to be counted with relative weight | or |. To be specific, all the spacelike plaquettes 
in the boundary, i.e. of the type Uij, get weight |. The timelike boundary plaquettes, Ui4, 
whose spacelike link lies in a face of the spatial cubic box of size [N — l^a^ also get weight 
|, while the timelike boundary plaquettes with spacelike link lying on one of the edges of 
this spatial box get a weight factor of |. This can be verified with a counting argument: 
There are 3N[N — 1)^ spacelike and 3N'^[N — 1) timelike plaquettes per time slice in 
the box with the boundary included. Of the spacelike plaquettes, 6(A^ — 1)^ lie in the 
boundary. Of the timelike plaquettes, 12(A^ — 1)^ lie in the boundary, 12(A^ — 1)(A^ — 2) of 
them with spacelike link in a face and the remaining 12(A^ — 1) with spacelike link on an 
edge of the [N — l^a^ cube. Taking into account the weight factors for these boundary 
plaquettes, we find the expected net number of 3(A^ — 1)"^ spacelike and 3(A^ — 1)"^ timelike 
plaquettes contributing to the energy, at a time-slice of the box. 

In the continuum, the monopole energy in a spatial box of volume [N — l^a^ is 
computed by subtracting the energy contained in the region outside this box from the 
infinite volume value (1) for the energy of the classical monopole. For large Nafi the 
outside contribution can be estimated by the integral of the asymptotic energy density 
l/g^r"^ over three-space, excluding the three-dimensional box. Denoting the length of the 
cubic box by 2R = [N — l)a and rewriting the integrand using 



1 dt 

—fexp[-r'^t], (31) 



„4 



t 



we get 



^out = ^(/ d'x--l (/3x-)=^/ ds[l - eie{s)] 




gm Jo 



, a = 1.2031 . (32) 
(/^ an 

Thus we obtain for the energy inside the box 

E,n = E- ^out = ^ f/" - ^ V (33) 

By comparing this result with the computed plaquette energy for various lattice sizes, at 
fixed a/i, and extrapolating to large /ii?, we can estimate the discretization error in the 
energy. For a/i = 7r/8, this discretization error is below the 1% level. Fven for a/i as large 
as 1.0 the discretization error is only 2%. 

Another check of discretization effects is provided by cooling (minimization of the 
action) of the latticized classical monopole, keeping the fields in the boundary fixed at 
their classical values, to find the corresponding local minimum of the lattice plaquette 
action. On an = 16 lattice, for example, the energy of an a/i = 1.0 monopole got 
lowered less than 0.1%. We conclude that the lattice version of the classical monopole 
is smooth enough on reasonable lattices such that discretization errors in the plaquette 
energy are small. 
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3.2 Boundary conditions on the lattice 

As discussed in Section 2.2, boundary conditions are imposed by fixing the gauge field 
components tangential to the boundary, cf. Eq. (18). On our lattice, this is a very natural 
procedure, see Fig. 5. We fix the spacelike links lying in the boundary of the cubic spatial 
volume and the timelike links originating at sites in the boundary, at each time slice. The 
boundary conditions in the time direction are taken periodic. 

The links lying in the boundary are fixed according to the appropriate link matrices 
(28,29), with exponential terms neglected. For the spatial links, this means K — > 0, as in 
(18). For the timelike links, we set H/r — > fi — l/r. Due to the loss of spherical symmetry, 
r varies in the boundary, so the result is different from (18), where /i — 1/r = /io was 
fixed. Actually, we shall set H/r — > /ii — 1/r, in order to be able to distinguish between 
the boundary parameter /ii and the induced monopole scale /i. We shall also study 
boundary conditions with H/r — > /io, i.e. with fixed /io- These may lead to approximately 
the same results if r may be considered effectively constant in the boundary, such that 
1^0 = 1^1 — 1/ -Reff- 

Summarizing, we shall distinguish two kinds of qualitatively different boundary con- 
ditions. We will call them BPS boundary conditions [K = 0, H/r = [ii — 1/r), 



k, tanet 



1 



r{x 



rr: ^akixi 



r[x 



(34) 



and HP boundary conditions [K = 0, H/r = /io) 



k, tang 

Al 



1 



r[x 

/io Xa 



(35) 



These names are inspired by the behaviour of H/r for the massless Bogomol'nyi-Prasad- 
Sommerfield solution (8,11) and the massive-Higgs version of the 't Hooft-Polyakov mono- 
pole (12), respectively. Also in this cubic volume the tangential boundary conditions are 
not incompatible with electric screening. 

If we adopt HP boundary conditions (35), the boundary links become 



2[X 



n{ X 



U^x) 



cos mix) ■ (7 sm 



\5 



1 + 



r^ -|- ya 



'■\/r'^ -|- 2ya -|- 



arctan 



wj x"^ ^ z 
r^ -|- ya 



^fx^ 



m ■ a , 



\5 



r^ -|- ya 



'■\/r'^ -|- 2ya -|- 



(and similarly for \J\ and 

aj-iQ . ^ _^ , (2/io 
cos I X ■ (7 sm . 



(36) 
(37) 

(38) 
(39) 

(40) 
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Here x{^) the angle between the vectors pointing from the origin to the beginning and 
end of the link. In the case of BPS boundary conditions (34), the spatial links are the 
same as in (36-39), whereas the timelike links are obtained by replacing fiQ ^ fii — l/r[x) 
in (40). 

For HP boundary conditions the timelike plaquettes in the boundary are exactly zero. 
For U24 for example, 

Ul{x, y, z) Uiix, y, z) U2{x, y, z) = Ui{x, y + 1, 2;) . (41) 

This agrees with the continuum result that the electric field is screened in the HP case, 
and hence zero if exponential effects are neglected. Similarly, one can show that for BPS 
boundary conditions the timelike plaquettes equal 1 — cos(a/2ri — a/2r2), independent of 
/io- Here ri^2 are the values of r at the sites between which the plaquette is located. Since 
the argument of the cosine is very small, (9(a^/i?^), this electric energy density is nearly 
equal to the continuum value oc Ijr'^. 

We have to assess the role of the shift symmetry (21) in our cubic lattice geometry. In 
the continuum it comes from the gauge transformation (20) and its twisted (antiperiodic) 
generalization. The normal (periodic) gauge transformation (20) can be taken over on 
the lattice without change. However, on the lattice we cannot use the twisted version 
because the link variables U^{x) are not invariant under the centre of the gauge group. 
The role of twisted gauge transformations is taken over by (normal) gauge transformations 
combined with the operation of multiplying each timelike link on a particular time slice 
by an element of the centre ( — /), which is a symmetry of the action and the integration 
measure. 

The transformation (20) applied to a smooth continuum configuration introduces a 
singularity in the spatial origin. On our lattice, the origin is located in the centre of an 
elementary lattice cube, so the only remnant of this singularity is that the link fields close 
to the origin get far away from the unit matrix. However, on the lattice we need not restrict 
ourselves to smooth gauge transformations. We can in fact restrict the transformations 
to the boundary, using the identity transformation in the interior. 

On the lattice, this gives rise to the same exact shift symmetry (21) as in the contin- 
uum, whether we adopt HP or BPS boundary conditions. This is guaranteed by the fact 
that = (j)2 = at the boundary, in the terminology of Fq. (14), so that (j) is invariant 
under the transformation, which acts on it as a rotation. (For the same reason, the abelian 
symmetry (19) of the boundary conditions also holds on our lattice.) The transformation 
of Aq at the boundary is given by 

- + <^ps), (42) 

f^o + -^ (HP), (43) 

while Ai does not enter into the boundary conditions. In both cases, the transformation 
induces a constant shift in /io or /ii, as in (21), even for the BPS case where the value of 
Aq depends on the boundary coordinates. 

Next, we consider the role of the anti-electric monopole. It is obtained from the 
monopole by mapping H/r — > —H/r^ at each value of x. In the case of HP boundary 
conditions, this means that /io is mapped into — /io in the boundary, so that we end 



/il 



r X 



1^0 
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up with HP boundary conditions for a different /io-value. Consequently, we envisage 
behaviour similar to that displayed in Fig. 4, the monopole having the same energy as the 
corresponding anti-electric monopole. For BPS boundary conditions, the mapping to the 
anti-electric monopole means /ii — 1 /r[x) — > — /ii -|- 1 /r[x) in the boundary. Because of the 
sign change in the 1 /r[x) term the anti-electric monopole does not obey the BPS boundary 
conditions (34). However, in sofar as BPS boundary conditions are approximated by HP 
boundary conditions, we expect that anti-electric monopoles can also be induced by BPS 
boundary conditions. This will be illustrated by the following computer simulations. 

We took an 8"^ X 12 lattice, and looked for local minima of the action for HP or BPS 
boundary conditions at given values of /io, /Ui- A random ('hot') starting configuration 
was created and its action was iteratively minimized ('cooled') subject to the boundary 
conditions. The actions of the resulting configurations were plotted as a function of /io and 
/ii. This was done for five different hot starting configurations, with identical boundary 
conditions, to allow different local minima to appear. The procedure was carried out for 
HP boundary conditions with /io ranging from to 2tt /T and for BPS boundary conditions 
with fii between and Att/T. 

Fig. 6 shows the resulting scatter plots of the energy (i.e., action per unit of time). 
The five points plotted at each value of /io, /<i are distributed among several local minima 
of the energy. We would like to identify them with the classical monopole solutions of 
the various symmetry-related scales discussed earlier. Thereto, we shifted the /io and fii 
values towards the appropriate values /ig and /i'^, according to (42,43), and replotted the 
data in Fig. 7. For HP boundary conditions (a), where fi^ and — /ig give the same energy, 
the points are plotted as a function of l/igl. The BPS data (6) are plotted against fj,[ itself. 

A striking feature in Figs. 66 and 76, for BPS boundary conditions, is the apparent 
refiection symmetry around fj,iT/2TT 0.5 (mod 1) and fj,[T/2TT 0.5, respectively. This 
suggests that the configurations on the descending branches of the curves in these figures 
can be interpreted to a good approximation as anti-electric monopoles. A quadratic 
fit demonstrates that the symmetry is almost perfect, see Fig. 76. The central value 
is fj,[T/2TT = 0.481, which means afj,[ = 0.252. We note that Figs. 7a and 6 almost 
coincide when afj,[ is shifted over 0.252, as is illustrated in Fig. 8. This confirms our 
expectation that, as far as the energy is concerned, HP boundary conditions are to a 
good approximation equivalent to BPS boundary conditions, with 

/"o = /"i — 1/-Reff (mod 27rn/r), (44) 

where a/i?eff = 0.252. Since i?eff i?, because it appears to be independent of /<oi, it 
follows from R/a = {N -l)/2 = 3.5 that 

i?eff = 1.13i?. (45) 

Figs. 6a and 7a, for HP boundary conditions, confirm the picture sketched in Sec- 
tion 2.2, cf. Fig. 4. Because of the HP boundary conditions together with the cubic 
geometry, we expect 1/R corrections to the asymptotic behaviour E-m oc /io, instead of 
only exponential corrections as in the cylindric case, cf. Fq. (24). In fact (omitting the 
primes on the /i's now), 

^„P, , ^Rps, . 47r, 1 , 47r / 1 1 \ 47r / 0.051\ 

El'{,o) - i^P„-(.0 = ^(.. - ^) = 7 (/^o + ^ - = ^ (.0 + , 

(46) 
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Figure 6: Energy, in units of Air/ag'^^ of local minima found on the 8"^ X 12 lattice as a 
function of [IqT I2tt, where T = 12, for HP boundary conditions (a), and as a function of 
fiiT/2TT for BPS boundary conditions (6). For each value of /io, /"i, five points are plotted, 
several of which coincide. 
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Figure 7: The energies of Fig. 6, now plotted as a function of [i'qT I2tt (a), or fj,[T/2TT (b), 
where obtained from /io, /"i through the appropriate shift (42,43). The HP case 

is depicted in a, BPS in b. For HP boundary conditions (a), the monopole (with /Xq) and 
the anti-electric monopole (with — /Xq) have both been plotted at positive /Xq, confirming 
that their energies are equal. In addition to the total energy (□), the magnetic (O) and 
electric (A) contributions are shown separately. The dashed line in 6 is a quadratic fit. 
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and the 0{1/R) correction to the energy formula (24) for the cylindrical case is small. It 
is also satisfying that the minimum [g'^ /Air) aE-m ~ 0.40 in Figs. 6-8 can be understood 
from (25) in terms of R^r- 1.627 a/i?eff = 0.41. 

We conclude that, with the appropriate modifications and the introduction of i?eff, the 
analysis in Sect. 2 can be taken over to the cubic lattice. 




Figure 8: The data of Figs. 7 a and 6, combined into one graph (total energy only), now as 
a function of a/ig (a few points have been left out). The BPS data (□) have been shifted 
over a distance l/i?eff along the horizontal axis. The HP points are plotted as diamonds 
(O) here. 



4 Monte Carlo simulations for the quantum mono- 
pole 

The quantum theory is described by the euclidean path integral 

Z[(3, {[/bnd}] = / DUint exp[-S{(3, {[/b„d}, {U-.^t})] • (47) 

Here {[/bnd} and {//int} are the sets of fixed boundary links and dynamical interior links, 
respectively. The integration symbol stands for a group invariant integration for each 
interior link U-mt- This path integral is a gauge-invariant functional of the link fields in 
the boundary. This is reminiscent of the gauge-invariant background field dependence of 
the path integral and the effective action in the background field method. 

We have seen that the monopole boundary conditions induce classical monopoles as 
local minima of the action in a (2i?)"^ X T lattice hypercube. Recall that R = a[N — l)/2. 
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T = aN4 on our x N4 lattice. The appropriate boundary conditions for the pure 
gauge theory are BPS, given in (34). It turned out that these are equivalent, to a good 
approximation, to HP boundary conditions (35). We may therefore restrict ourselves to 
HP boundary conditions in the dynamical simulations to follow. This is fortunate since 
they would be more appropriate, if our assumption is correct that the classical dyon turns 
into a purely magnetic monopole in the quantum theory. 

As already mentioned in Sect. 2, for the Monte Carlo computations of the quantum 
monopole the global minimum of the action is important. For the 8"^ X 12 example in 
Fig. 6 this corresponds to the envelope of the lowest data points. It follows that we may 
restrict the values of /io to < /io < tt/T. 

Figs. 6,7 show rather large finite box effects for these values of /io- Classically we 
can eliminate such effects by using the exact BPS values for the boundary conditions 
(including 'exponential effects'), as in Sect. 2.1, but in the quantum theory we do not 
know the precise form of the monopole. This is why we use the universal long distance 
form of the monopole fields in formulating the BPS and HP boundary conditions. To 
reduce the finite box effects we shall use the largest value /io = vr/T in the numerical 
simulations. Note that even then fioR is only 0.92 for the 8"^ X 12 lattice in Fig. 6. 

Having chosen HP boundary conditions for the Monte Carlo simulations, characterized 
by the boundary parameter /io, we note that the /i-scale of the monopole depends on the 
nature of the quantum monopole. If it stays BPS-like, then /i = /io + l/i?eff, as in the 
classical theory. If, on the other hand, it becomes HP-like in the quantum theory, then 
its /i-scale is /i = /io. So we have 



a/i = a/io + 1/(1.13 * (iV - l)/2), /ii? = (^/2)(iV - l)/iV4 + 1/1-13 (BPS), (49) 



according to the quantum monopole turning out HP-like or BPS-like. In the classical 
regime, at weak coupling, we know the monopole behaves BPS-like. At stronger coupling, 
we will have to keep both options open for our analysis of the results of the simulations. 

The finite box effects can be decreased by increasing /io-R, i.e. increasing the ratio 
N/N4. However, we will have to deal with fioR values that are not very large, for the 
following reasons. Firstly, in the quantum case we should not take R too large if we per- 
form nonperturbative computations. The reason is that nonperturbative effects generate 
a finite correlation length, as given by the glueball mass or the string tension, implying 
that boundary conditions become irrelevant for large volumes. To stay in the semiclassical 
regime we have to use a finite volume, with R\/g smaller than some number z of order 1. 
In the monopole description of confinement this restriction can be interpreted by saying 
that we would like to study one monopole, which may include the quantum fiuctuations of 
other small scale objects, such as instantons or monopoles, but which should exclude non- 
perturbative effects of other large scale monopoles, as these destroy the monopole under 
study. Secondly, we want to study monopoles with inverse scales up to ^/a 1. Since 
^i'"^ \/a < fiQ^ ^/o^ < z/fioR, we have a conflict between large fioR and getting |J~^^/o^ 1. 

We recall that the value of the lattice distance a, e.g. in units of the string tension -^0", 
0-^0", is controlled by the bare coupling constant (3 = 4/(/^. On a given X N4 lattice, 
decreasing (3 from inflnity means increasing R\/o and fj,~^y^ from zero onwards, keeping 
fiR constant. In this way we can access sizes up to R\/g of order 1, while it depends on 
the value of /ii? if we can also reach fj,~^y^ 1. 
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Here we shall use a symmetric lattice, N = A^4, as a compromise between getting 
sufficiently large fj,~^y^ and not too small fiR. The resulting values of fiR are rather 
small, cf. (48,49). Consequently our Monte Carlo simulations deal mainly with the interior 
of the monopole. For HP, fiR is even smaller than for BPS. On the other hand HP-like 
monopoles are expected to be smaller, at given /i, than BPS-like monopoles, so their 
interior may still fit reasonably well in the box. 

4.1 Strategy for calculating the monopole mass 

We write the monopole mass as 

M = Mi„ + Mout , (50) 

where Min is the contribution to the mass which will be computed in the Monte Carlo 
simulation and Mout comes from the infinite region outside. As discussed above, the latter 
cannot be dealt with nonperturbatively, it makes sense only in an analytic semiclassical 
approximation. Here we shall take it into account as follows. 

In the pure SU(2) gauge theory the classical monopole is BPS, but in the quantum 
case we shall explore the possibility that it turns HP-like and analyse our results both 
for BPS and for HP behavior. For a classical BPS monopole the electric and magnetic 
fields decay as 1/r^ at large distances. For an HP monopole the magnetic field decays as 
1/r^ while the electric field is exponentially small. In the box the outside contribution 
to the energy is 1/aR (BPS) or l/2aR (HP), in units of 47r/(/^, where a = 1.2031 is the 
correction factor introduced in Fq. (32). In the quantum theory the outside magnetic field 
is expected to remain that of a magnetic monopole of charge —1, whereas the expected 
electric field corresponds to electric charge 1 in the BPS case, and to charge in the HP 
case. 

Thus we write for the outside contribution to the monopole mass 

Mo„t(/i) = J^^, . CoutifiR), Cout = (BPS), (51) 



1 



(HP) . (52) 



2a/ii? 

Presumably it would not matter very much if we use the renormalized coupling in other 
reasonable schemes such as MS or Pauli-Villars, because the corresponding A scale pa- 
rameters are very close, and furthermore (51,52) is supposed to be a correction to the total 
mass. On the other hand, our calculation of Min is scheme independent, and a change in 
the definition of the coupling is compensated by a change in C\n- 
We define the monopole mass Min for the interior region by 

Min = 4ln%i^. (53) 

Here Znion(/5; o/io) is the partition function (47) for a monopole background subject to 
boundary conditions specified by the parameter [Iq and Zvac(/3) is the analogous definition 
for vacuum boundary conditions, i.e. all = 0. Fq. (53) defines the monopole mass as 
the free energy of the monopole, with the vacuum free energy subtracted. In this way, the 
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exponentiated monopole mass expresses the effective Boitzmann weight of the monopole 
in the path integral. In our simulations we will also monitor the internal energy AE", 

= ^{{SUon - (5')vac), (54) 

and its density function of r = 

We assume the mass (53) to be ultraviolet finite. Renormalization properties of func- 
tional integrals depending on boundary conditions are studied in [21], but the results 
there are not directly applicable to our case. Reassuring to us are also the following ob- 
servations. The are finite in perturbation theory in the sense that the divergent parts 
of the wavefunction and coupling constant renormalizations cancel. (We use the conven- 
tion in which 'contains' the bare coupling constant g. We have of course employed 
the same convention for the boundary values, e.g. /io does not require renormalization). 
Furthermore the path integral Z^^n depends continuously on the fields in the boundary. 
These boundary fields can be sent to zero in a continuous fashion, leading Z^^n towards 
Zvac and the monopole mass to zero. (Also in the continuum the boundary fields can 
be deformed to zero continuously.) Note that in the limit /3 — > oo the monopole mass 
assumes its classical value. 

The mass (53) can be expressed in terms of expectation values of observables, making 
it accessible to Monte Carlo computation. Differentiating Eq. (53) we find (recall S oc 
(3 = W) 

This leads to 

Mi„(/3) = f^mP) , (56) 

where the integration constant corresponds to Min = at /3 = 0. This is clear from S oc (3 
and the normalized measure for the integration over the interior links. 
We use the trapezoid rule in computing 

N ^ 

aM,^{fi) = aJ2 2 (A " A-i)(A^(A) + A^(A-i)) , (57) 

8 = 1 

for a set of values /3o = < /3i < • • • < (^n-i < f^N = /5 such that (57) is a good 
approximation to (56). This requires two Monte Carlo simulations at each value of 
one for monopole and one for vacuum boundary conditions. 



4.2 The internal energy 

We shall present the results of a simulation on an 8"* lattice. Recall that this means our 
spatial box has volume (2i?)"^ = (7a)"^ and the length of the periodic time direction is 
T = Sa. For convenience we collect here the numerical values of some of the quantities 
introduced earlier. We have Res/a = 1.13 i?/a = 3.96. For HP, a/i = tt/S = 0.39 and 
liR = 1.37. For BPS, a/i = tt/S + a/R^s = 0.65 and jiR = 2.26. 

The simulations were carried out as follows. We fixed the boundary conditions ap- 
propriately and took a hot (random) start in the interior of the lattice. We started at 
the lowest /3- value, then stepwise increased /3 to /32, /^s, etc. up to /3jv, and subsequently 
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we decreased (3 again, down to (3i. This allowed us to check for possible "hysteresis" 
effects, due to unsufficient thermalisation. As a matter of fact, we carried out two such 
cycles. The second one was done to obtain data at intermediate /3-values where more 
accuracy appeared necessary after the first run. See the specifications detailed in Table 
f. At each /3- value we started with 400 thermalisation sweeps and subsequently 4000 
measurements of the action were done, separated by 2 sweeps. The statistical errors on 
AE were determined by pointwise subtracting the 2 sets of 4000 numbers (using the same 
pseudo-random numbers for monopole and vacuum boundary conditions) and binning the 
4000 differences. The errors on aM-m were determined similarly. 

Our SU(2) program used the site labelling and "gather and scatter" method described 
in Ref. [22], and additional logical vectors, used as masking vectors throughout the pro- 
gram, controlled the special behaviour of the links in the boundary, see Fig. 5. The 
pseudo-heatbath method [23, 24] combined with uj = 2 overrelaxation [25] in the ratio 1:1 
was used for updating. 

The data for the internal energy AE are given in Table 1. The data for increasing 
and decreasing (3 are statistically compatible, in other words no "hysteresis" effects are 
observed. Therefore, we will present graphs for the l/cr^-weighted averages only. 

We note that aAEjfi contains a constant contribution 



leading to a term cfi in aM\n. For the 8"* lattice c = 0.0380. This contribution is purely 
magnetic and comes entirely from the part of the fixed plaquettes in the boundary that 
is attributed to the energy inside the box, cf. the discussion after Fq. (30). In the limit 
/3 = the links become uncorrelated so that the interior does not contribute to AE/(3. 
For the presentation of the data for AE and Min, however, it will be convenient to ignore 
the contribution (58). Therefore we will suppress it until we come to the calculation of 
the physical monopole mass and the function C in Section 4.3. 

In Fig. 9 AE is displayed as a function of (3. For /3- values up to 1.8, AE differs little 
from zero. In this region the correlation length is very small and the boundary conditions 
do not cause an observable effect in the internal energy, within errors. 

The steep increase at /3 2.1 signals that the boundary conditions manage to induce 
a monopole in the box. As a consequence of the increasing correlation length the interior 
of the box is then no longer decorrelated from the boundary. It is possible that the 
appearance of the monopole coincides with the box transition, the analogue of the finite 
temperature transition for a symmetric hypercube, which is at (3c ~ 2.33, for a 7"* lattice 
with periodic boundary conditions [26]. The shift from (3 = 2.33 to 2.1 could be due 
to the difference in boundary conditions. However, it is also possible that the box/finite 
temperature transition in Z^^n and Zyac is not observable in AE. This can be studied in 
large boxes ^ A^4, where one would expect the appearance of the monopole at a /3mon 
much larger than (3c. 

For larger /3, the internal energy slowly decreases towards its asymptotic value, deter- 
mined by cooling 20 random configurations satisfying the monopole boundary conditions, 
cf. Section 3.2. This asymptotic value is slightly smaller than the value following from 
Fig. 7a, for a/io = 7r/8 ~ 0.39 {ii'qT I2tt = 0.75 in the figure), because the contribution of 
the plaquettes lying in the boundary is omitted. The physically most interesting region 
extends from f3 = 2.3 onwards, where one expects scaling behaviour. 



c = aAE I f3\p=o , 
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0.041 ±0.042 


0.085±0.081 


-0.025±0.055 


0.017±0.039 


-0.008±0.084 


1.90 


0.128±0.044 


0.083±0.050 


0.211±0.082 


0.100±0.040 


0.056±0.038 


0.156±0.070 


1.95* 


0.104±0.048 


0.048 ±0.040 


0.152±0.081 


0.136±0.038 


0.070±0.038 


0.206±0.070 


2.00 


0.158±0.050 


0.098±0.050 


0.257±0.097 


0.117±0.044 


0.071 ±0.040 


0.188±0.080 


2.05* 


0.159±0.051 


0.123±0.041 


0.282±0.088 


0.121±0.040 


0.071 ±0.038 


0.191±0.072 


2.10* 


0.252±0.036 


0.150±0.034 


0.403 ±0.062 


0.355 ±0.044 


0.230 ±0.045 


0.585±0.085 


2.15* 


0.430 ±0.039 


0.292±0.032 


0.722±0.067 


0.394±0.055 


0.288 ±0.045 


0.683±0.090 


2.20 


0.457 ±0.030 


0.334±0.031 


0.791 ±0.052 


0.468 ±0.030 


0.319±0.030 


0.786±0.053 


2.25* 


0.494 ±0.030 


0.357±0.026 


0.851 ±0.047 


0.428 ±0.030 


0.290±0.029 


0.718±0.051 


2.30* 


0.406 ±0.026 


0.275±0.026 


0.681 ±0.043 


0.435 ±0.030 


0.282±0.027 


0.717±0.049 


2.35* 


0.461 ±0.026 


0.309±0.024 


0.770 ±0.040 


0.461 ±0.022 


0.305±0.027 


0.766 ±0.040 


2.40 


0.406 ±0.030 


0.248 ±0.025 


0.653 ±0.045 


0.410±0.023 


0.239±0.027 


0.650 ±0.042 


2.50* 


0.417±0.018 


0.278±0.021 


0.695±0.030 


0.359±0.019 


0.227±0.021 


0.586±0.030 


2.60 


0.371±0.018 


0.218±0.020 


0.589±0.030 


0.393±0.023 


0.190±0.020 


0.583±0.037 


2.70* 


0.366±0.016 


0.199±0.021 


0.564±0.030 


0.385±0.022 


0.185±0.020 


0.570±0.032 


2.80 


0.378±0.015 


0.199±0.016 


0.577 ±0.024 


0.354±0.015 


0.176±0.016 


0.530±0.023 


2.90* 


0.367±0.014 


0.176±0.016 


0.542 ±0.022 


0.381±0.015 


0.185±0.018 


0.566±0.025 


3.00 


0.367±0.016 


0.175±0.016 


0.542 ±0.024 


0.373±0.015 


0.186±0.015 


0.559±0.023 


3.20* 


0.364±0.013 


0.153±0.015 


0.517±0.021 


0.360±0.012 


0.167±0.013 


0.527±0.019 


3.50 


0.364±0.011 


0.151±0.013 


0.516±0.018 


0.361±0.010 


0.141±0.012 


0.503±0.017 


4.00 


0.368±0.011 


0.150±0.011 


0.518±0.014 


0.352±0.010 


0.156±0.010 


0.507±0.015 


CXD 


0.3841 


0.1116 


0.4957 





Table 1: Monte Carlo data for AE", in units of Air/ag'^^ from two cycles of measurements. 
The magnetic and electric contributions separately as well as the total are given. In the 
first cycle (3 was gradually increased from 0.2 to 4.0 and subsequently decreased back to 
0.2. The second cycle, indicated by the * on the (3 values in the Table, ran from (3 = 1.95 
to 3.2 and back. Each AE value in the table is based on 4000 configurations. The zero 
entries for (3 = refiect the fact that the contribution AE / [Air / ag'^)\p=o = c/tt = 0.0121 
of Eq. (58) is not included in this table. The (3 = oo entry is the classical value determined 
by minimization. 
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Figure 9: The internal energy AE in units of Air / ag'^ as a function of (3 = 4/(/^. Shown are 
the magnetic (O) and electric (A) components and the total (□). For (3 < 2.1 only the 
total is plotted for clarity. The dotted horizontal lines indicate the asymptotic (/3 — > oo) 
values, determined by minimization, cf. Table 1. 

The magnetic and electric components of the internal energy are seen to follow essen- 
tially the same behaviour as a function of (3. The fact that the asymptotic, classical, value 
of the electric energy is much smaller than that of the magnetic energy is due to the small 
value of fioR (cf. Fig. 7a at [i'qT I2tt = 0.75, which corresponds to the same value of fioR). 
The magnetic and electric internal energies do not differ as much in the intermediate (3 
region as in the limit /3 — > oo. As a matter of fact, the magnetic energy approaches its 
limiting value very slowly, and from below, whereas the electric contribution approaches 
its asymptotic value from above. 

It is also interesting to look at the spatial density profile of the internal energy. It 
is shown in Fig. 10 for several values of (3. For clarity we present it without error bars 
and with data points connected by straight lines to guide the eye. The energy density 
is given in units of the inverse lattice spacing a. However, since a/io is fixed at 7r/8, 
and R/a = 3.5, this is equivalent to units of either of the monopole scales [i = [Iq ov 
= /io + l/i?eff- The energy density at r is computed by averaging the energies of all 
elementary cubes whose centre is at a distance r from the centre of the box. The energy 
in such a cube is proportional to the sum of the plaquettes forming it. The resulting r- 
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dependent distribution is averaged over time slices and Monte Carlo configurations. The 
figure shows clearly that a lump of energy builds up when the monopole is induced in the 
box for (3 > 2.2. With increasing /3, the energy density decreases towards the classical 
/3 — > oo curve, in accordance with the observation made for the total internal energy, 
Fig. 9. The wiggles in the classical curve give an indication of the size of discretization 
effects and of boundary effects due to the cubic geometry. 




Figure 10: Fxpectation values for the internal energy density in units of 47r/a'*(/^, as a 
function of r/a. 

The behaviour of the energy density as a function of (3 does not provide conclusive 
information as to whether the quantum monopole becomes HP-like or stays BPS like. 
(The data for electric and magnetic contributions separately are also inconclusive.) Note, 
for comparison, that the finite temperature electric screening mass ^2/3 gn/aN^ is fairly 
small for couplings gj^ ^ 10. 

4.3 The monopole mass and the function C{gji) 

Now we turn to the computation of the monopole mass. The results for Min in units of 
the lattice spacing are shown in Fig. 11. These data have been obtained by integrating 
the data for the internal energy AE" of Table 1 and Fig. 9, using Fq. (57). The constant 
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contribution c(3 from the boundary plaquettes (58) is still omitted. At each given (3- 
value, the integration was carried out for all the 4000 configurations separately, allowing 
an estimate of the statistical errors on the computed masses by binning. Because of the 
integration procedure, the masses at neighbouring values of (3 are strongly correlated. 




Figure 11: The dimensionless monopole mass aM-m as a function of /3, for a/io = 7r/8. 
Apart from the total mass (□), the magnetic (O) and electric (A) contributions are 
shown separately. For (3 < 2.1 only the total mass is plotted for clarity. 

In the limit /3 — > oo, the mass curve approaches a straight line with slope given by 
aAE[(3 — > oo) = 0.4957 x 47r/(/^ = 1.5573/3, see Fig. 9. This line is approached from 
below because AE is a decreasing function at large (3. A fit of the mass data in the large-/? 
region, including the c(3 contribution, gives 

aM-JfJ) ^ 1.5953 (fJ - 1.415) - 0.97 ^ , (59) 

P 

with prescribed slope 1.5573 + 0.0380 = 1.5953. 

We now use the mass data to compute the function C in the mass formula Fq. (2). 
For the interior of the box we write 

47r Q/jj/ 

aMin{j3; a/io) = 2/a / a CinisU^R/ f^j) ^ (60) 
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where we do the analysis for the HP and BPS values of /i given in (48,49). In order to 
extract C-m from the mass data, using this formula, we need to compute Ar/ fi for both 
choices of /i and subsequently gj^[AR/ fi) as a function of (3. Then the total C is obtained 
by taking into account the correction (52) or (51) for the outside region. 

We use lattice string tension data to determine the lattice distance a in terms of a 
physical scale. We write 



— = a/i^^— . (61) 
Ar aJa Ar 



The quantity \/o j Ar is an input parameter, since the Monte Carlo simulations have not 
yet determined it unambiguously. The data on the string tension a^/a were collected 
from various sources [27-29]. We interpolated some of them towards neighbouring /3- 
values. The values we used, tabulated in Table 2, come from Polyakov-loop correlations 
and include the finite-size correction 7r/3iv^ to a [30-32]. Possible systematic errors on a 
are inherited by our function C . We have analysed our data for three different values of 
y/a I Ar: 1.7, 2.0 and 2.2. The former is consistent with both the value given in Ref. [33] 
and the value obtained by combining the fi = 2.85 result a^/a = 0.0636 of Ref. [29] with 
the lower bound 0.038 on gAr given there. The value 2.0 for / Ar follows from Ref. 
[34], and 2.2 is the variational estimate we obtained in Ref. [1]. 

The running coupling gj^[AR/ fi) is determined from the Ar/ fi values thus obtained by 
means of the two-loop /3-function for this coupling, 

- -f3{9R) = --^ = bogl + Wgt, , (62) 

with coeflicients bo = ll/(247r2), bi = (102/121) 6^ for SU(2). The exact solution at this 
order in gj^ is given by 

t = (4 + 1) ° (54) • 

We find g\{ARl [i) by inverting Eq. (63) numerically, instead of resorting to an approxi- 
mate analytic inversion. The use of the two-loop formula becomes questionable at distance 
scales of the order l/ y/a, at the 'border of the semiclassical region', when becomes 
large. We keep in mind, therefore, that this may cause a systematic error of perhaps 20% 
in our HP-results at /3 = 2.3. 

This analysis leads to the values for C given in Table 2, where we have included 
the contribution cfi = 0.0380/3 from the plaquettes in the boundary (58). The quoted 
numbers take into account C^^ = 0.302 or C^^^ = 0.368, from (52,51) and (48,49). The 
/3 = 00 entries were obtained from Cin = A£'/(7r/3/i), [i — > 00, giving = (0.49577r -|- 
0.0380)/(0.393^) = 1.293, C^^ = 1.596, and CfJ'^ = (0.393/0.646) = 0.787, C^^^ = 
1.155. These C- values are useful for reference, although at /3 = 00 the classical pure SU(2) 
monopole is BPS with C = 1. 

Fig. 12 shows C as a function of gj^[AR/ fi). Since //bps > /"hp the points in the BPS 
analysis are found at weaker running coupling than those in the HP analysis. The dotted 
horizontal lines are the /3 — > 00 values 1.596 and 1.155 described above. The error bars 
correspond to the statistical errors on the mass data only. The data points in each set 
should lie on a smooth curve because they are strongly correlated, as in the case of the 
mass data of Fig. 11, and the error bars represent the interval within which these curves 
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/3 




aMin 


^jkn = 1.7 


y^/AR = 2.0 


y^/AR = 2.2 


9r 


c 


9r 


c 


9r 


C 


2.20 


0.467 


0.609 ±0.063 


15.46 


1.95±0.20 


11.42 


1.52±0.15 


9.99 


1.37±0.13 


2.25 


0.415 


0.735 ±0.063 


12.26 


1.92±0.16 


9.70 


1.58±0.12 


8.69 


1.45±0.11 


2.30 


0.368 


0.856 ±0.064 


10.25 


1.90±0.13 


8.46 


1.62±0.11 


7.70 


1.50±0.10 


2.35 


0.313 


0.973 ±0.064 


8.47 


1.82±0.11 


7.26 


1.60±0.09 


6.71 


1.50±0.09 


2.40 


0.266 


1.084±0.066 


7.25 


1.76±0.10 


6.37 


1.58±0.08 


5.95 


1.50±0.08 


2.50 


0.186 


1.288 ±0.064 


5.58 


1.65±0.07 


5.06 


1.53±0.07 


4.81 


1.46±0.06 


2.60 


0.132 


1.482 ±0.066 


4.60 


1.59±0.06 


4.25 


1.49 ±0.06 


4.07 


1.44±0.05 


2.70 


0.101 


1.670 ±0.067 


4.06 


1.59±0.05 


3.79 


1.51±0.05 


3.65 


1.46±0.05 


2.80 


0.0742 


1.850 ±0.067 


3.58 


1.57±0.05 


3.37 


1.49 ±0.05 


3.26 


1.45 ±0.04 


2.90 


0.0545 


2.027±0.067 


3.21 


1.55±0.04 


3.04 


1.48 ±0.04 


2.95 


1.45 ±0.04 


3.00 


0.0400 


2.205 ±0.067 


2.91 


1.53±0.04 


2.77 


1.48 ±0.04 


2.70 


1.45 ±0.04 


CO 





CO 





1.596 





1.596 





1.596 




/3 


a^fa 


aMin 


v^/Ar = 1.7 


y^/AR = 2.0 


y^/AR = 2.2 


9r 


C 


9r 


C 


9r 


C 


2.20 


0.467 


0.609 ±0.063 


7.70 


0.95±0.06 


6.70 


0.87 ±0.05 


6.24 


0.84±0.05 


2.25 


0.415 


0.735 ±0.063 


6.94 


1.00±0.05 


6.14 


0.92 ±0.05 


5.75 


0.89±0.04 


2.30 


0.368 


0.856 ±0.064 


6.33 


1.04±0.05 


5.66 


0.97 ±0.04 


5.33 


0.93±0.04 


2.35 


0.313 


0.973 ±0.064 


5.66 


1.05±0.04 


5.13 


0.98 ±0.04 


4.87 


0.95±0.04 


2.40 


0.266 


1.084±0.066 


5.13 


1.05±0.04 


4.70 


1.00 ±0.04 


4.48 


0.97±0.04 


2.50 


0.186 


1.288 ±0.064 


4.27 


1.05±0.03 


3.97 


1.00±0.03 


3.82 


0.97±0.03 


2.60 


0.132 


1.482 ±0.066 


3.69 


1.04±0.03 


3.47 


1.00±0.03 


3.35 


0.98±0.03 


2.70 


0.101 


1.670 ±0.067 


3.34 


1.06±0.03 


3.16 


1.02±0.03 


3.06 


1.00±0.03 


2.80 


0.0742 


1.850 ±0.067 


3.02 


1.06±0.02 


2.87 


1.02±0.02 


2.79 


1.00±0.02 


2.90 


0.0545 


2.027±0.067 


2.75 


1.06±0.02 


2.63 


1.03±0.02 


2.57 


1.01±0.02 


3.00 


0.0400 


2.205 ±0.067 


2.53 


1.06±0.02 


2.43 


1.03±0.02 


2.37 


1.01±0.02 


CO 





CO 





1.155 





1.155 





1.155 



Table 2: Computed values for gj^ and C, for the HP analysis (top) and the BPS analysis 
(bottom). The analyses have been done for the three quoted values of y/a j Kr, as described 
in the text. The values for the string tension used for the calculation, obtained from 
Refs. [27-29], are given as well. 
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Figure 12: C as a function of gj^[AR/ fi). The different sets of points are for / Ar = 1.7 
(O), 2.0 (□) and 2.2 (A), for both the HP analysis (upper set of curves) and the BPS 
analysis (lower set). The eleven points for each value of y^/ Ar and /i are based on the 
mass data for (3 = 2.2, 2.25, 2.3, 2.35, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, from right to left. 

as a whole can move. The small deviations of the data points themselves from the smooth 
behaviour of the dashed lines in Fig. 12, drawn to guide the eye, are due to errors in the 
string tension data used. 

The solid curves in Fig. 12 give a rough indication of the approach to the classical 
limit. They have been obtained from the fit (59). Since there are no string tension 
data available for large /3, we used the 2-loop approximation to the /3-function for the 
bare lattice coupling here to compute a A/, and hence Ar/ fi for both values of /i, using 
Ar/Al = 20.78 [10]. It is well known that this is a bad approximation for intermediate 
/3, so these curves should be viewed with caution. 

In the small-(/|j region the monopole will be BPS-like. The deviation of the classical 
C- value of 1.155 from the true value C = 1 of the classical BPS monopole is primarily 
an 'exponential correction' due to the small value or fiR in our simulation. To a lesser 
extent, it is also caused by our use of HP boundary conditions for simulating the BPS 
monopole (cf. the slight difference between HP and BPS in Fig. 8). Thus, the small-(/|j 
points in the HP analysis in Fig. 12 are misleading. At larger values of gj^ we have two 
scenarios. If the monopole remains BPS-like, one of the lower curves applies, depending 
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on the value of y/a j Kr. If the monopole becomes HP-like, on the other hand, the graph 
of C will lie in between the two sets of curves in the plot, starting off at C = 1 and rising 
to one of the upper curves at intermediate or large g\. 

In the HP scenario, the conclusion we draw is that the mass formula (2) holds over a 
wide range of g\ values. Using only values of fi in the scaling regime fi > 2.3, and taking 
a BPS monopole in the /3 — > oo limit as discussed above, we see that the function C is 
slowly increasing as a function of In the BPS scenario, we had no conjecture about 
the behaviour of C . (Recall that we introduced this function because of the analogy with 
the mass formula for the 't Hooft-Polyakov monopole.) Thus, there is no reason to expect 
it to be an increasing function of g\. Nevertheless, even in this case C varies little, albeit 
over a smaller range of running couplings. 

5 Conclusion 

We have performed a numerical investigation of the magnetic monopole in pure SU(2) 
gauge theory. The classical monopoles could be caught as local minima on a lattice with 
prescribed boundary conditions. Gauge symmetries limited the scales [i of the monopoles 
corresponding to global minima, relevant for the Monte Carlo simulations in the quantum 
theory. The preferred scale was related to the finite extent of the system in the temper- 
ature direction. Since the boundary conditions induce the quantum monopole only for 
finite volume, we were led to studying the interior of the monopole in a finite temperature 
situation. The exploratory numerical data for the interior of the monopole show how 
the energy density (Fig. 10) and mass (Fig. 11) behave as a function of the bare gauge 
coupling. 

In view of our earlier work [1] we analyzed these data with two scenarios: the quantum 
monopole could remain BPS-like or become HP-like. Our data for the energy density were 
insufficient to decide between the two alternatives. Assuming that the monopole behaves 
HP-like at intermediate and large values of the coupling, we found the function C((/|j) 
in (2) to be slowly increasing, for renormalized coupling values ranging from g\ ?> to 
8 (C ?^ 1.6), or 3 to 10 (1.6< C<2.0), depending on the value of ^/a j Kr (cf. Fig. 12). 
These coupling values correspond to scales [i varying from about lO-^o" to ^/(T, which 
include the [i of order y/a needed in Ref. [1]. 

As a variant of the HP interpretation we can take the simulated monopole at face 
value. It is an object with magnetic charge, and a free energy given by the mass formula. 
The shape of the energy density suggests that these monopoles keep their identity when 
packing them at a minimal distance h = fj,~^ . 

Fven in the BPS scenario the mass formula appears to have reasonable validity. The 
corresponding function C was found to decrease only slightly between = and 5.5 or 
7. 

Hence, in both scenarios the scale dependence of the monopole mass is essentially 
given by the scale dependence of the renormalized coupling. This supports our assumption 
about the monopole mass mentioned in the introduction. 

These results can of course be improved by going to larger lattices. In particular, 
controlled matching with the classical limit C = 1 for gR — > 0, which was not possible 
with our small spatial volume, may be possible using asymmetric X N4 lattices with 
N ^ N4. This may also allow a more detailed investigation of the screening properties of 
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the electric charge. 
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